Penetration depth study of LaOs4Sbi2: Multiband s-wave superconductivity 
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We measured the magnetic penetration depth A(r) in single crystals of LaOs4Sbi2 (rc=0.74 K) 
down to 85 mK using a tunnel diode oscillator technique. The observed low-temperature exponential 
dependence indicates a s-wave gap. Fitting the low temperature data to BCS s-wave expression 
gives the zero temperature gap value A(0) = (1.34 ± 0.07)fc_BT'c which is significantly smaller than 
the BCS value of 1.76kBTc. In addition, the normalized superfluid density p{T) shows an unusually 
long suppresion near Tc, and are best fit by a two-band s-wave model. 



The skutterudite PrOs4Sbi2 has been receiving much 
attention due to the discovery of its unconventional heavy 
fermion superconductivity, and the possible role played 
by quadrupolar fluctuations in its pairing mechanism i^r— 
An important step in clarifying its superconductivity is 
to determine the symmetry of the superconducting gap. 
Conflicting results had been reported. On one hand, 
scanning tunneling microscopy (STM), muon spin rota- 
tion (/xSR), and nuclear quadrupolar relaxation (NQR) 
experiments point to an isotropic superconducting energy 
gapi^^-? On the other hand, magnetic penetration depth, 
angle-dependent thermal conductivity, and neutron scat- 
tering measurements revealed point nodes on the super- 
conducting gap4^— Angle-dependent thermal conductiv- 
ity data revealed two distinct superconducting phases of 
twofold and fourfold symmetries, both with point nodes 
on the superconducting gap.^ Recently, multiband super- 
conductivity was proposed as a reconciliation for these 
conflicting results. According to field-dependent thermal 
conductivity data, PrOs4Sbi2 is a two-band supercon- 
ductor, with nodes on one band and an isotropic gap on 
the other. ^"'^"^^ A similar two-gap structure was observed 
in Li2Pt3B, the nodes there arising from the interference 
of singlet and triplet pairingJ^ 

LaOs4Sbi2 (superconducting transition temperature 
rc=0.74 K) and PrOs4Sbi2 (Tc=1.85 K) are iso- 
structural superconductors. The substitution of La by 
Pr introduces the 4/ electrons, while largely preserv- 
ing the crystal structurei^"— and Fermi surface topol- 
ogyJ^ The 4/ energy levels of the Pr^+ ion are split 
by the crystal electric field (CEF), resulting in a non- 
magnetic Pi singlet ground state, a 0.7-meV Ps first ex- 
cited state, and 11-meV P4 second excited state (in Oh 
cubic symmetry).— The higher Tc of PrOs4Sbi2 com- 
pared to LaOs4Sbi2 is attributed to the type of scat- 
tering between the conduction electrons and these low- 



lying CEF-split 4/ levelsi^"— — the inelastic scattering 
of the conduction electron by the Pi— ?>P5 transition has 
a strong quadrupole matrix element that enhances pair 
formation and consequently Tc- This interaction, known 
as aspherical Coulomb scattering, overcomes the mag- 
netic pair breaking effect due to the s-J exchange scat- 
tering of the conduction electron by the Pi^'r4 transi- 
tion, which has a strong dipole matrix element that sup- 
presses Tc- The / electrons thus play an important role 
in skuttcrudite/heavy- fermion superconductivity via the 
CEF-split /-electron energy levels. 

Since it has been shown that PrOs4Sbi2 is a two-band 
superconductor with different pairing symmetries on each 
band, we want to find out if LaOs4Sbi2 is a multiband 
superconductor as well, and if so, whether the symme- 
tries of these multiple gaps are different. Such an inves- 
tigation of LaOs4Sbi2 will allow us to determine whether 
the /-electrons are necessary for the formation of multi- 
band superconductivity. In general, multiband super- 
conductivity arises when there are more than one con- 
duction band and the pair couplings within each band 
are of different magnitudes 1^. Though filled skutteru- 
dites generically contain two hole-like conduction bands, 
observations of multiband superconductivity were only 
reported for PrOs4Sbi2 and PrRu4Sbi2 — both contain- 
ing /-electrons.^'' These /-electrons might exert differ- 
ent degrees of influence on the pair coupling within each 
band, hence giving rise to multiband superconductivity. 
In this context, it is tempting to attribute multiband su- 
perconductivity in Pr-based skutterudites to the presence 
of /-electrons. 

There have been evidences that LaOs4Sbi2 is a weak- 
coupling s-wave superconductor. Specific heat mea- 
surement showed a discontinuity jump AC at Tc with 
AC/-fTc = 1A6^ close to the BCS value of 1.43. Sb- 
NQR measurement showed a coherence peak at Tc and 



exponential temperature dependence at low tempera- 
turesii The data can be fitted by an isotropic gapli^ 
with A(0) « f.GfcsTc or an anisotropic gap A{6) = 
<5 + (A - S)sme with A = f .TSfesTc and S ^ f .2f /csTc^ 
However, to the best of our knowledge, multiband super- 
conductivity has never been reported for this material. 

Magnetic penetration depth has been shown to be 
a valuable tool for probing multiband superconductiv- 
ity, for instance, in the well-known multiband super- 
conductor MgB2.^^ At low temperatures, the deviation 
AA(T)=A(T) — A(0) is sensitive to the low energy exci- 
tations of quasiparticles across the superconducting gap. 
Consequently, the symmetry of the energy gap can be 
determined from the temperature-dependence of AA(T). 
In addition, fits to the superfiuid density over the entire 
temperature range of measurement allows one to detect 
the existence of multiple gaps as well as their temper- 
ature dependences. In this Paper, we report the mea- 
surement of magnetic penetration depth A(T) in single 
crystals of LaOs4Sbi2 down to 85 mK. We found that at 
low temperatures, AA(r) exhibits an exponential tem- 
perature dependence consistent with a s-wave gap sym- 
metry. Fitting the data to BCS s-wave expression gives a 
minimum gap value that is significantly smaller than the 
BCS value, suggestive of multiband superconductivity. 
The superfiuid density data over the entire temperature 
range are best fitted by a two-band s-wave model. 

Details of sample growth and characterization are de- 
scribed in Ref. IS. Measurements were performed utiliz- 
ing a 21-MHz self-resonant tunnel diode oscillator™ with 
a noise level of 2 parts in 10^ and low drift. The sam- 
ple was positioned in the center of a induction coil which 
forms part of the LC resonant tank circuit. The AA(r) is 
directly proportional to the change in resonant frequency 
A/(T) of the oscillator, with the proportionality factor 
G dependent on sample and coil geometries. For a square 
sample of side 2w, thickness 2d, demagnetization N, and 
volume V, G is known to vary as G oc i?3_D(l — N)/V, 
where Rzd = w/[2{l + {l + 2d/wf)a.rcia.n{w/2d)-2d/w\ 
is the effective sample dimensioui^ For our sample 2w k, 
0.59 mm and 2d w 0.13 mm. The value of G was deter- 
mined for a high-purity Al single crystal by fitting the Al 
data to the extreme non-local expression, and then ad- 
justing for relative sample dimensions. The magnitude of 
the ac field inside the induction coil was estimated to be 
less than 40 mOe, and the cryostat was surrounded by a 
bilayer Mu-metal shield that reduced the dc field to less 
than 1 mOe. The single crystal sapphire rod which held 
the sample with GE varnish was attached to the mixing 
chamber of an Oxford Kelvinox-25 dilution refrigerator 
to provide cooling. During the experiment, the probing 
ac field of the inductor coil was directed along the c-axis 
of the sample and the in-plane penetration depth Xab was 
measured. Since LaOs4Sbi2 has a cubic crystal structure, 
we will omit the distinction between axes for the rest of 
this paper. 

As A(0) w 4700 A in LaOs4Sbi2r^ it probes a signif- 
icant depth into the sample and is therefore less sensi- 



2 




0.08 0.12 0.16 0.20 0.24 
T(K) 



FIG. 1. Temperature dependence of the penetration depth 
AA(r) in sample 1 at low temperatures. The solid line is a fit 
to Eq.([T]). Inset: AA(r) over the entire temperature range. 

five to surface quality, giving result representative of the 
bulk. At low temperatures (T < Tc/S), it is well estab- 
lished that A(T) of an isotropic s-wave superconductor 
will asymptomatically approach an exponential behavior 
given by 

where A(0) and A(0) are the zero-temperature values of A 
and the superconducting gap, respectively. For a weakly 
anisotropic gap or multiple gaps, Eq. ([T]) will still follow 
but now with A(0) being replaced by the minimum gap 
value in the system, and A(0) by an effective value which 
depends on the details of gap anisotropy. 

We measured the AA(r) in two single crystals of 
LaOs4Sbi2. An exponential temperature dependence 
characteristic of a s-wave energy gap is observed in the 
two samples over the temperature range 85 mK-0.26 K. 
Figure [T] shows the low-temperature variation of AA(r) 
in sample 1. The inset in Fig. [1] shows AA(T) of sample 
1 over the entire temperature range. The AA(r) starts 
to fall at 0.76 K and reach the transition mid-point at 
0.68 K, consistent with the reported value of rc=0.74 K 
in other experiments 

The low-temperature data of sample 1 was fitted to 
Eq. HI) up to 0.25 K (~Tc/3), as shown by the sohd 
line in Fig. [T] The zero-temperature gap value A(0) ob- 
tained from the fitting was found to be (0.99±0.05) K or 
(1.34±0.07) fcsTc (assuming 7;=0.74 K), where the er- 
ror corresponds to fitting the data up to T = 0.23 K and 
T = 0.26 K respectively. For sample 2, fitting the low- 
temperature data to Eq. (P) gives A(0) = (1.06±0.5) K or 
A(0) = (1.43 ± 0.07)fcBrc, slightly larger than the value 
obtained in sample 1. Wc also fitted the low-temperature 
data to a power law AA(r) ^ T" and obtained n—A, ex- 
cluding the possibility of line nodes in the energy gap 
which gives AA(r) T or AA(T) in the presence 

of impurities. As sample 1 has a sharper transition near 
Tc, we will focus on it for the rest of this paper. The 
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value A(0) = 1.34fcsT'c is significantly smaller than the 
weak-coupling BCS value of UGfe^Tc. This suggests the 
possibility of anisotropic gap or multiple gaps (multiband 
superconductivity) . 

The scenario of multiple gaps is more likely, since in 
addition, we also observed an unusually long suppression 
of the normalized superfluid density p{T) — A^(0)/A^(T) 
near Tc which could not be due to gap anisotropy (Fig.[2|). 
On the contrary, the presence of multiple gaps often re- 
sults in suppression of p(T) near due to the presence of 
the smaller superconducting gap .^^i'^° This feature should 
be intrinsic to the material since it has been observed in 
another two samples and could not have been caused by 
some magnetic ordering, as /iSR experiment observed no 
spontaneous magnetism even down to 20 mK.^^ 

We propose a two-band model to fit our data. In this 
model, the Hamiltonian is given by 



ik(7 ikk' 
kk' 

where i=l,2 represents the first and second bands. 
ci^ka and C2,fecr are the corresponding electron opera- 
tor. Vi^kk' = Vi7i,fc7i,fc' and V2,kk' = V'272,fe72,fc' arc 
the reduced pair coupling for the two bands. V^^kk' = 
V37i,fe72,fe' is the intcrband pair coupling. This model 
has been used to describe the two-band superconductor 
MgB2.3i 

Taking the BCS mean-field approximation, the inter- 
action between the two band is decoupled as 



H - ^ f-i,kctkcr(^i,ka + 



ika 



-fci ■ 



A* E7fcC, _fe^c,,fct (3) 



where Ai = Y.k (^i7i,fc(ci,-fe;ci^fe^) + "1^372, fc(c2,-fc4,C2,fet)) 
and A2 = Y.k (^272,fe(c2,-HC2,fct) + ^37i,fe(ci,-fciCi,A;|)). 
We propose the tight-binding band dispersion 

U.k = — 2ii (cos kx + cos ky + cos kz) 

i kx cos kz -\- cos ky cos k. 



-|-4i ■ (cos kx cos ky -\- cos / 
-Ml- 



(4) 



Here ti{i — 1,2) are the hopping constants, t'^ are the 
next nearest-neighbor hoppings, and /Xi are the chemi- 
cal potentials of the two bands. The same tight-binding 
band dispersion with a slight modification due to an- 
tiferromagnetism has been used to describe multiband 
superconductivity in the electron-doped cuprate super- 
conductor Pr2_3;Cea;Cu04i^ 

We next solve, self-consistently, a set of equa- 
tions of Ai,A2,/ii,M2 obtained by SF/5x = 
{x = Ai, A2, Ml, M2), where F = -^InZ = 
— ■glne^'^^ is the free energy, for given model param- 
eters ti,t'^,Vi, ¥3,121, 122- The superfiuid density along a 
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FIG. 2. Experimental data (O) and theoretical fit (solid line) 
for the normalized superfluid density p{T). Insert: tempera- 
ture dependence of the two energy gaps from the theoretical 
fit. 

certain direction, such as 2:-direction, is given by 



EE 

1=1,2 fe 



' dkl 



dkx 



where f (Ei^k) is the Fermi function. In our fit- 
ting, we assume two bands to have s-wave gaps, 
namely, 7i,/c=72.fe=l- The choice of gap sym- 
metry is supported by the low-temperature ex- 
ponential behavior in AA (this work) and Sb- 
NQR experiment.'' For LaOs4Sbi2, the model pa- 
rameters used are {niti,t[,Vi,n2t2,t'2,V2,V3) — 
(0.08, t2/4.4, -0.25t2, 0.084, 0.25, ^2, -b.3i2, 0.103i2, 0.00 
08t2) . Here the hopping integral <2 = 1 is taken as the 
unit of energy in the whole fitting, and is related to 
critical temperature as Tc « 0.049t2. 

From the fitting, we obtained Ai(0) = 1.69kBTc, 
A2(0) = l.SOksTc, and 712/711 = 3.13. The calculated 
value of A2(0) — 1.30fc_BTc agrees well with the value of 
the zero-temperature gap value A(0) = (1.34±0.07)A;bTc 
obtained from low-temperature data of AA(T). The two 
values agree, since at low temperatures, the smaller en- 
ergy gap will be the effective gap for quasiparticle exci- 
tation. The calculated value of Ai(0) = l.GdksTc is also 
consistent with the weak-coupling value of l.lQksTc^^ 
Figure [2] shows that the experimental data of p(T), in 
particular the long suppression near Tc, can be well fit- 
ted by our two-band model. The long suppression is effec- 
tively caused by the smaller gap A2 vanishing at higher 
temperatures due to interband coupling with the bigger 
gap Ai (Inset of Fig. [5]). 

Our present finding reveals a two-band s-wave super- 
conductivity in LaOs4Sbi2 . Given the close Fermi topolo- 
gies of PrOs4Sbi2 and LaOs4Sbi2, the present result 
shows that multiband superconductivity of PrOs4Sbi2 
persist even when the /-electrons are removed. This ex- 
cludes the /-electrons from being the origin of multiband 
superconductivity in PrOs4Sbi2. Though the /-electrons 
are not responsible for the multiband superconductivity, 
they may however influence the behavior of each band. 
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Our present finding of s-wave symmetry in botli bands 
suggests tliat the two bands evolve differently in the pres- 
ence of /-electrons, namely into a nodal band and a 
fully-gapped band in PrOs4Sbi2. This unconventional 
effect of /-electrons in PrOs4Sbi2 may be related to the 
weakly split CEF levels,— since it was found that the 
nodal feature disappears once PrOs4Sbi2 is doped with 
32,33 therefore suggest that theoretical models for 
PrOs4Sbi2 should take the multiband superconductivity 
of LaOs4Sbi2 as a starting point, and develop the sub- 
sequent modifications in superconducting properties due 
to the adding of /-electrons. 

The magnetic penetration depth data of PrOs4Sbi2 
and PrRu4Sbi2 had been measured by one of us with 
the same experimental technique employed here^i^!^^ 
We tried to fit the data of PrOs4Sbi2 with a similar 
two-band s-wave model, but was not successful — this is 
consistent with the unconventional (point node) nature 
of one of its gaps. Also, we were able to fit the PrRu4Sbi2 
data using just one superconducting gap, showing that 
the multiband superconductivity of PrRu4Sbi2, as shown 
in thermal conductivity data, does not significantly affect 
the magnetic penetration depth data. 



In summary, we report measurements of the magnetic 
penetration depth A(T) in single crystals of LaOs4Sbi2 
down to 85 mK. We find that A(T) exhibits an exponen- 
tial temperature dependence at low temperatures with 
the zero-temperature gap value A(0) ~ (1.34±0.07)fcsrc- 
Our results show that LaOs4Sbi2 is a two-band s-wave 
superconductor. Given the close Fermi topology of 
PrOs4Sbi2 and LaOs4Sbi2, this suggests that the /- 
electrons are not the origin of multiband superconduc- 
tivity in these two materials. 
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